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Abstract. Yudin's lower bound 1211 for the spherical designs is generalized 
to the cubature formulas on the projective spaces over a field K C {IR, C, H} 
and thus to isometric embeddings 1™ K — > Z™. K with p S 2N. For large p and in 
some other situations this is essentially better than those known before. 
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1. Introduction 

In the theory of spherical and projective designs some important lower bounds 
were obtained [H [3] by maximization of the functional 

£>(/) = feKlt 

c o[f\ 

where Ki is the set of nonncgativc on (—1, 1) nonzero polynomials /, degf < I, and 

C0[/]= j f(t)LJ a ,p(t)dt, LO a fi(t) = (\-t) a (l+tf. (1.1) 

Here the numbers I € N and a, f3 > — 1 depend on the design. 

Obviously sup {£>(/) : / g K{\ = sup {/(l) : / g K h c [f] = 1}. The solution to 
the latter linear programming problem is classical, the extremal polynomial / max is 
unique and can be expressed in terms of the Jacobi polynomials, see |19j . Section 
7.7.1. For the designs of cardinality n this yields 

n > Ta,/3/ma X (l) = T a £ max {£>(/) : / g Ki} , (1.2) 

where 

T a .fi = / U a fi{t) dt. (1.3) 



We denote by L^' 0) (-l, 1) the space of complex- valued measurable functions / 
on (—1,1) such that 

2 _ 



EE J \f(t)\'cJ a , (t)dt<OO. 

The corresponding Jacobi polynomials Pk(t) constitute an orthogonal basis in 
L^"' \ so that 

oo 

f(t) = J2^[f]Pk(t), feL^l (1.4) 

fc=0 

where 

Ck[f}= [ f(t)P k (t)uj a fi(t)dt, ^. = l/||P fc || 2 . (1.5) 
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(|1 .4|) The Jacobi-Fourier series ()1.4|) converges to / in (— 1, 1)- The coefficient 

Co[/] in (|1 .5|) coincides with that of since -Po(i) = 1, according to the usual 

standardization 



degP fc = fc, P fe (l) = 




(k = 0,1,2,...)- 



For the same reason vq = l/r a)/ g. 

The linear programming bound (| 1 . 2[) can be extended to the set Kiji, I' > I, of 
the polynomials / ^ 0, deg/ < Z', such that f(t) > for |i| < 1 and Cfe[/] < for 
I + 1 < k < V . In this way Boyvalenkov and Nikova [31 HI [5] obtained a series of new 
concrete lower bounds for the projective designs. For the spherical designs Yudin 
[2"T] considered the limit case V = oo. Its class Ki t00 consists of all nonnegativc 
nonzero continuous functions f(t), \t\ < 1, such that Ck[f] < for all k > I + 1. A 
suitable choice of a function / € Ki tOQ yields a lower bound asymptotically better 
than classical one that comes from (|1.2ll . 

In the present paper we generalize Yudin's result on the projective designs and 
even the cubature formulas on the projective spaces KP m_1 , where K C {K, C, H}. 
(Recall that H is the standard notation for the quaternion field.) The extension of 
the linear programming bound from the projective designs to the general cubature 
formulas is technically simple but important since the latter are equivalent to the 
isometric embeddings Z™k ~~ * ^kj P £ 2N. (See [14] and references therein.) Note 
that with the standard inner product (x, y) the space Z™k ^ s Euclidean, its unit 
sphere is S = S Sm ~ l , where S = 5{K) such that 6(R) = 1, 5(C) = 2, 5(B) = 4 . 

From now on we assume rn > 2, p G 2N, and denote by $k(to,p) the space 
of complex- valued functions 4>(x), x £ S, satisfying the following conditions, see 

Hang. 

(1) (f> — ip\S, where ip is a homogeneous polynomial of degree p ("p- forms") on 
the space K m = R Sm ; 

(2) ip is invariant in the sense that 

ip(wa) = ip(w) (w£K m ,aeK, \a\ = 1). 

A fortiori, <p(xa) — 4>(x) that allows us to naturally transfer (j) to the projective 
space KP m_1 . However, we will consider cj> on S which is equivalent but more 
elementary. In this setting a projective cubature formula of index p on S is 

/n 
4>dcr = ^2<f)(xi)pi, 0€$ K (m,p), (1.6) 

i— 1 

where a is the normalized Lcbcsgue measure on S, the nodes x% € S are projectivcly 
distinct and the weights pi arc positive. (Note that pi = 1 automatically by the 
restriction of (x,x) p ^ 2 to S.) In the case of equal pi the set {x^}™ is nothing but a 
projective p/2-design, c.f. [9j. 

2. Basic theory 
First of all, we have the decomposition 

p/2 

$ K (m,p) = ^ Harm K (m, 2k) (2.1) 

fc=0 
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where the space HarrriK(m, 2k) consists of restrictions to 5* of the invariant harmonic 
2/c-forms. Regarding to the inner product 

(i>i,i>2) = J i>ii>2 da 

the decomposition (|2.ip is orthogonal. 

For any orthonormal basis {4>ki}i={ 2k of Harms (m, 2k) the addition formula 

d m ,2k 

^ (f>ks(x)(f> ks {y) = b m ^ k P k {xy) (x,y € S) (2.2) 

s=l 

S(m -I) -2 n 8-2 . . 

m,k = T at pv k P k (l), a = , = — , (2.3) 



holds with 



and 

xy = 2\(x,y)\ 2 -l, (2.4) 

see [SI [T31 LH]. Later on we operate only with a, /? given by (|2.3[) . 

Now let X be a finite nonempty subset of S, and let A(X) be its angle set, i.e. 

A(X) = {xy : x,y G X,x £ y). 

The addition formula easily implies the following 

Lemma 2.1. Let the series 

oo 

^a k P k {t) 

k=Q 

converge to a function f(t) for every t € A(X) and for t = 1. Then 

OO 2k 

J2 f{xy)\{x)\{y)=Y j akb- m ) k £ 

x,y£X k=0 s=l 

where A is an arbitrary function X — > C. 



fcs (a;)A(a:) 



(2.5) 



With A(x) = 1 formula (|2.5[) plays a fundamental role in the design theory 
j6[[T0l[T3j. In the context of cubature formulas we need (|2.5|) with arbitrary X(x) > 
0, A(x) = 1, c.f.|T4|, §5. Also note that, in contrast to those which are quoted 
above, now we have to apply (|2.5p to the non-polynomial functions / € Ki tOC . It is 
possible because of 

Lemma 2.2. TTie J acobi- Fourier series of any function f G -?0.oo converges to f(t) 
for all t 6 [-1, 1]. 

Proof. Since /(£) is continuous, its Jacobi-Fourier series at t = 1 is summable to 
/(l) by a Cesaro method, see [T5], Theorem 9.13. Therefore, it is summable to 
/(l) by the Abel method, see [7], Theorem 5.5. Hence, this series converges to f(l) 
since c k f < for k > I + 1. It remains to refer to Theorem 7.32.1 from [19] which 
states that 

max|P fe (i)| =P fe (l) (2.6) 
l*l<i 

if max(a, /3) > -1/2. The latter is fulfilled because of ([23]) and m > 2. □ 
Corollary 2.3. Formula (|2.5j) is true for every f £ Ki tOQ with a k = v k c k [f], k > 0. 
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Remark 2.4. In [21j the absolute convergence of the corresponding series is men- 
tioned without proof. The proof of Lemma 12.21 shows that in our situation the 
convergence is absolute and uniform. 

Now we can prove the following linear programming bound. 

Proposition 2.5. The inequality 

n > r a . 3 sup {£>(/) : / G K p/2j00 } (2.7) 

holds for any projective cubature formula of shape (|1.6|) . 

Proof. We have 

n 

^2<f>(xi)ei=Q, (f>e Harm K (m, 2k), l<k<p/2. (2.8) 

»=i 

Applying Corollary |2 . 31 and Lemma [2~T1 to / G K p / 2oo , X = {xt}™ and X(x) = p(x), 
x G X, we obtain 

i£X x,y£X 



Indeed, on the left side of ()2.5[) all summands are > 0. On the right side the 
summands arc < for k > p/2 + 1 and vanish for 1 < k < p/2 by (|2.8p . It 
remains to recall that p(x) = 1, therefore, J2p 2 ( x ) — n ~ X \ on the other hand, 
a o&mo = c o[/]/^a,/3 since 6 m , = 1. □ 



Remark 2.6. The inequality (|2.7|) implies 

n > r Qi/3 sup{£(/) : / G K p/2J ,}, I' > p/2, 

since K p / 2 ,i> C K p/2t00 . (For Z' = Z we set = Jfj.) 

Since with any given m,pa projective cubature formula exists (or, equivalently, 
there exists an isometric embedding IV? K — > Z" K ), we have 

Corollary 2.7. sup{L>(/) : / G 

The supremum in question is unknown but a "good" test function can be con- 
structed using the "convolution" 



g(xu)h(uy) da(u) (x,y €z S) (2-9) 
of two suitable functions g(t) and h(t), — 1 < t < 1, c.f. [2Tj . 

Lemma 2.8. For any e G 1, 1) the function u 1— ► e(su), u € S, belongs to 

Li(S, a) for every x € S and 

[ e(xu)da(u) = — [ e(t)u a p(t) dt (2.10) 

J T a jj J_ l 

Proof. This follows by calculation in spherical coordinates. □ 
Corollary 2.9. With g,h G L 2 ^\— 1, 1) £/ie integral (|2.9p exists for all x, y G 5. 
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Since any ordered pair x' , y' € S with x'y' = xy can be obtained from x, y by an 
isometry of Z™ K , the integral (|2.9|) depends on xy only. Thus, we have a function 
(g * h)(t), -l<t< 1, such that 

(g*h)(xy) = J g(xu)h(uy)da(u). (2.11) 

In particular, for x = y (|2.1ip yields 

(g * ft)(l) = / g{xu)h{xu) da(u) = — f g{t)h(t)u a .p{t) dt, (2.12) 

J T a,P J-l 

by (|2.10j) . Moreover, applying the Schwartz inequality to (|2.11| and using (|2.10|) 
again we obtain 

sup|(ff*ft)(i)| < — \\g\\-\\h\\. (2.13) 

By this inequality and bilinearity, the convolution g * h determines a continuous 
mapping (L^ ,,3) ) 2 -> L M . 

Lemma 2.10. ITit/i 6 the function (g * h)(t) is continuous, and the 

series 

oo 

(<? * h){t) = Y / (^k/b m ,k)c k [g}c k [h}P k (t) 

k=0 

converges uniformly. 



Proof. Let 



Then 



since 



N N 

g N = ^2 v i c ] Id] p j i h N (t ) = ^2 Vk Ck W Pk ■ 

j=a k=o 

N 

g N *h N = y^^l/b m ,k)ck[g]ck[h}Pk 

k=0 



Pj * Pk = b-^PkSjk 

by the addition formula. Since g^ — ► g and Hn — * ft (-ZV — > 00) in L^'^, we obtain 
<?iV * ftAT g * h uniformly Thus, the limit function is continuous. □ 

Corollary 2.11. c k [g * ft] = ffcC fc [g]c fe [Zi]/& mife = c k [g]c k [h]/T a ^P k {\). 

3. A FUNCTION /j G A'/.oo 

Recall that all roots of every P k are simple and lie on (— 1, 1). The roots of the 
derivative P' k alternate them, so they are also simple and lie on (—1,1). Now we 
introduce a function // by setting 

/,=.,., *>-{*<•>-*<«• ;« m 

where r = Z + 1 and £ is the largest root of P' r . We have to verify that // € Ki t00 . 

By Lemma \2. 101 f; is continuous. The inequality fi >0 follows from (|2. 1 1[) since 
ft > and g > 0. The former is obvious, the latter is true since <?(£) = 0, g(l) > 
and g'(t) ^ for f < t < 1. Moreover, / ; (1) > by (|2~T21) . thus, /j ^ 0. It remains 
to prove that c k [fi] < for k > r. In [5T] a rather complicated vector analysis on 



(i 



YU. I. LYUBICH 



R m was used at this point. We manage without a generalization of this technique 
to C m and H m by dealing with the corresponding Jacobi polynomials. 
Our starting point is the differential equation 

^ i = (u a+1 ,p+ 1 Piy + i(i + X)u) ai pPi = 0, i>0, (3.2) 

where A = a + p+l, see [19], formula(4.2.1). Note that A > by (pOj) , From (f3~2|) 
it follows that 



= J (P r A k - P k A r )dt 



(k-r)(k + r + A) 



y Lj a ,pP r p k dt + J {u a+h0+1 (p r Pi - p k p;.y dt = 



{k-r){k + r + X) / LJ a ,pP r P k dt - {uj a+U3+1 P r P' k ){i) 



since w Q+ i^+i(l) = 0, P/(£) = 0- For fc^rwe obtain 



(fc - r)(fc + r + A) ' 
This formula extends to r = since Po(t) = 1, so Pq(£) = 0- Thus, 



Pk^a,0 dt = 



k(k + X) 
and then 

/ PrPkVc.p dt= fc(fc *? P '® : ; / PkOJ a ,pdt. 
Je (k-r)(k + r+l) J £ 



As a result, 



r l f 1 p M r(r + X)P r (Q f 1 r(r + A)P r (Q 



and, by Corollary 12. Ill 



[k - r)(fc + r + X)P k {l)T a ^ 

Since Pfc(l) > 0, formula (|3.3[) yields signcfc[/;] = sign(P r (£)), k > r. But 
signP r (£) = —1 since £ lies in between two largest roots of P r (t) and P r (l) > 0. 
Thus, c k [fi) < for k > r. In addition, c r [fi] = since c r [h] = 0. The latter follows 
from (|3.2j) with t = r by integration over [£, 1]. 

In conclusion we note that £ in (|3.1|) is actually the largest root of P ; (t), 
see [19], formula (4.21.7). 

4. Main Theorem 

Now we are in position to prove the following 

Theorem 4.1. The number n of nodes of every projective cubature formula of 
index p on S 6 ™ 1-1 satisfies the inequality 

r( a + 2)r(/3 + i) /]\'(™-i>/» 

" r(o + /9 + 2)F(-/3,a+l,a + 2,e) \e J ' K ' } 
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where F is the hypergeometric function, the numbers a and [3 are given by (|2.3|) . 
£ = (1 — £)/2, £ is the largest root of the Jacobi polynomial P^J^'^ +1 ' (t) . 

Proof. Using f p /2(t) as a test function in (|2.7|) we get 

T a ,f3fp/2(1) 

n - — \f 1 ■ 

c 0Up/2\ 

By (|2~T2|) and (f3TT]) we have 



T a,/3/ P /2(l) = y ghLj a .fjdt = J guj a ,pdt = co[g]. 

On the other hand, co[/ p /2] = co[<?] c o[^]/ r a-/3 by Corollary 12. Ill Hence, 

n > °^ = ' \ ' (4 2) 

coW /^(l-tjo^ + t)/'* 

Now we substitute t = 1 — 2s into the numerator and t = 1 — 2e,s into thtc denom- 
inator. This yields (|4.1|) since 



a + 1, a + 2, e) = (a + 1) / s Q (l - es)' 3 <fs, 

Jo 

(c.f. [J, formula (15.3.1)) and 

Jo r(a + /J + 2) 

(Also note that a + 1 = (<5m - (5) /2 by (|23)) .) □ 

Remark 4.2. By substitution t = 2s 2 — 1 in both integrals (|4.2p we obtain 



n >JoV y 7(1 + 072. 

In particular, for K = M. we have a = (m — 3)/2, /3 = —1/2, see (|2.3[) . Hence, 

^(l-s 2 )^- 3 )/ 2 ^ 

n > ^ " (K = R), (4.3 

" f^(l- s 2 )( m - 3 )/ 2 ds 

where rj is the largest root of the polynomial 

p«--l)AV2) (2s 2 _ 1} = CQnst . p ((m-l)/2,( m -l)/2) (s)/s5 (4 4) 

or, equivalently, of the Gegenbauer polynomial C^i(s) (see [19], formulas (4.1.5) 
and (4.7.1).). In the case of antipodal spherical (p+ l)-design the lower bound (|4.3p 
turns into (3) of [5T] up to the additional factor 2 in the latter. Note that the factor 
2 is just the degree of the natural mapping S™ -1 — ► IRP m_1 . 

Remark 4.3. By ()2.3|) we have a = m — 2, j3 = for K = C, and a = 2m — 3, /3 = 1 
for K = H. Accordingly, (|4.1|) reduces to 

n > f-") (K = C), (4.5) 
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and to 



n > 



1 



(K = H). 



(2m- 1) - (2m-2)e \e 
hi the real case the hypergeometric function in (|4. 1 [) is not a polynomial of e 



(4.6) 



Now we denote by N^(m,p) the minimal number n of nodes in the cubature 
formula (|1.6[) or, cquivalcntly, the minimal n such that there is isometric embedding 
/ 



2:K 



/™. K . In this notation Theorem 14. II states that 



N K (m,p) > 



r(a + 2)r(/3 + 1) 



<5(m-l)/2 



(4.7) 



T(a + + 2)F(-P, a + 1, a + 2, e) 

We will compare this result to the linear programming bound (|1.2|) with ^ = p/2. 
An explicit form of the latter is 

N K (m,p) > A K (m, g), q = p/2, 



where 



while 



A K (m, g) 



V m-1 /' 

/•m+[g/2]-l\ /m+[( 9 +l)/2]-l\ 

V m-i A m-1 /' 

1 (2m+[q/2]-2^ /2TO+[(g+l)/2]-l^ 



k. 2m-l 



2m-2 



2m-2 



(K = K); 
(K = C); 
(K = H), 



Ak(to,p) < A K (m,p), 
(See [14] and the references therein.) 

5. Asymptotic analysis 
From and ijO]) it follows that 



(4.8) 



(4.9) 



(4.10) 



N K (m,p) 



> 



P 



5(m-l) 



Ak("z) ' 



P 



where 



2 m - 1 (m - 1)!, 
A K (m) = <( 2 4 (™- 1 )(m- l)! 2 , 

28(™- 1 )(2m- l)!(2m-2)!, 



(5.1) 



(5.2) 



or, in an unified form. 



, / x r(Jm/2)r(J(m-l)/2 + l) 25(m _ 1} Tja + (3 + 2)r(a + 2) 25(m _ 1) 

AkM = rp 2 = fcsTT) 2 

(5.3) 

As to (|4.7p . £ is the only parameter depending on p. (Of course, e also depends 
on m.) By definition, e = (1 — £)/2 = sin 2 (#/2) where = arccos£. This is 



the smallest root of the polynomial P 



(a+1,/3+1) 



p/2 



(cos#). By Theorem 8.1.2 from 



we have 9 ~ 2_7 Q ._|_i ) i/j!? where j Q +i,i is the smallest positive root of the Bessel's 
function J a+ i(z). Therefore, e ~ j 2 +1 j/p 2 , and (|4. 7[) yields 



N K (m,p) 



T(a + 2)T(P+ 1) p 5 ^- 1 ) 



oo, 



T(a + /3 + 2) 

v A ' Ja+1,1 

since e — > 0, F(-, •, -,0) = I. This estimate is better than (|5.ip because of 



(5.4) 
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Proposition 5.1. The inequality 

.a (m _i) r(a + 2)T(^+l) 

holds for all m > 2, except for the case m = 2, 5 = 1, w/ien (|5.5p changes for an 
equality. 

Proof. By (|5.3p the inequality (|5.5[) is equivalent to 

We set a + 1 = v, so that <5(m — 1) = 2u, and (|5.6p takes the form 

j£i<r(i/+l) 2 .16". (5.7) 

The number f is positive integer or half-integer. ^ > 1/2, and — 1/2 if and only 
if m = 2, 5 = 1. In this case j„,i = 7r since Ji/2( z ) is proportional to sm z/y/z. On 
the other hand, J?(3/2) 2 • 16 1 / 2 = 7r as well. Thus, (|5.7[) changes for an equality. 

Now let f > 1. By the inequality j u<1 < y / 2(v + l)(v + 3) (see [20], Section 15.3) 
it suffices to prove that 

(i/ + 1)"(i/ + 3)" <I> + 1) 2 -8". (5.8) 
By Stirling's lower bound the inequality (|5.8p follows from 

A fortiori, (|5.8[) follows from 

2tw f 4t^I > e 4 . 



But the latter is indeed true if v > vq where Vq is a unique root of the equation 
27ri/(8/e 2 ) iy = e 4 . It is easy to see that v < 6, so (|5.8[) is valid for ^ > 6. For v < 6, 
i.e. z/ = l,3/2,2,...,5,ll/2, the inequality (|5.8|) can be checked numerically. □ 



The inequalities (|5.ip and (|5.4p can be rewritten as 

limMp- 6( - m -^N K (m,p) > l/A K (m) (5.9) 

and 

,. . , -S(m-l) / N ^ r(a + 2)T(J) + 1) .-S( m -1) 

hminfp 1 >N K (m,p) > 3 a +i,i (5-10) 

p^oo 1 (a + p + 2) ' 

respectively. By Proposition 15.11 the quotient kk(to) of the lower bounds (|5.9p 
and (|5.10p is less than 1, except for the case m = 2, 5 = 1. Moreover, kk(to) 
exponentially decays asm-> oo, c.f. [21| for the spherical designs. 

Proposition 5.2. Asymptotically, 

1 (e\ s ( m - 1 ) 

«K(m) ~ — — - , m^oo. (5-11) 



7T<5m V4 

Proof. In the same notation as before we have 

i 2 " 

KK(m) = 



r(z/ + 1) 2 • 16" 
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Using the relation j v .\ = v + 0{y/v), v 
Stirling's asymptotic formula we obtain 



(see 



Section 15.83) and 



2ttv 
5m 



L(£Y 



k k (to) 

that is equivalent to (|5. 1 1(1 since 2v • 
Remark 5.3. From (|4.10[) the asymptotic upper bound 

limsnp p- 5 ^ m ^N K {m,p) < 2 s( - m -^ / X K (m) 



□ 



(5.12) 



follows. We see that there is an exponential gap between (|5.12[) and (|5. 10|) as 
77i — > oo. Indeed, the quotient of these bounds is 

1 /e\ 5 (™- 1 ) 
,2J 

2 



2^ m ^K K (m) 



00. 



(5.13) 



TV 5m 
6. The case m 

In this case we discuss the real, complex and quaternion situation separately. 



6.1. K = R. Then the inequalities (|4.7|) and (J478J) are both the equalities, so they 
coincide. Indeed, Ng.(2,p) = p/2 + 1, according to (T6J [18], and, on the other 
hand, Ar(2, q) = q + 1 = p/2 + 1 by (|4.9p . Furthermore, in the real case (|4.7p 
is equivalent to (|4.3|) . For m = 2 this yields N&(2,p) > 7r/2arccos7/ = p/2 + 1. 
Indeed, in this context r\ is the largest root of the Gcgenbauer polynomial Ci+i( s ) 
= sin(p + 2)9/sm8 where 9 = arccoss. 



1 



6.2. K = C. By (|49|) 

N c (2,p) > 

c.f.[Tl]. On the other hand, our bound (|4.5|) for 77 

iV c (2,j7)> 



2 is 



(6.1) 



(6.2) 



(t) and ]C[ means the smallest ingeger > C, 
coincides with (|6.1|) for p < 16, 



where ^ p is the largest root of P p j 2 
( £ 1. A numerical evaluation shows that 
but exceeds it for 18 < p < 90. Moreover, the difference Ac(p) between the lower 
bounds (|6.2[) and (|6.ip is nondecreasing in this range, as we see from the table 



p 


< 16 


18 - 24 


26,28 


30,32 


34, 36 


38 


40 


42, 44 


46 


48 


50 


52 


A c (p) 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 



P 


54 


56 


58 


60 


62 


64 


66 


68 


70 


72 


74 


76 


78 


80 


82 


84 


A c (p) 


12 


13 


14 


15 


17 


18 


19 


20 


22 


23 


25 


26 


28 


29 


31 


32 



P 


86 


88 


90 


Ac(p) 


34 


36 


38 



The table also shows that the "derivative " Ag(jj) = Ac(p) — Ac(p — 2) is 
nondecreasing (rather slowly). 



LOWER BOUNDS FOR PROJECTIVE DESIGNS, CUBATURE FORMULAS AND RELATED ISOMETRIC EMBEDDINGS 



6.3. K=H. We have 



Ak(2,,)>!(^ 



[(P + 2)/2] 
2 



from (gSj) and (|4~9j) . but (|4~6j) yields 

JVh(2,p) > 

,(2,2) 



(2 + 77 P )(l-77 P ) s 



(6.3) 



(6.4) 



where ry p is the largest root of (t) ■ 

Comparing (|6.4[) to (|6.3[) one can see a small advantage of (|6.3[) when 4 < p < 20, 
p 7^ 18. Namely, for the difference An(p) between the lower bounds (|6.4[) and (|6.3[) 
we have 



p 


2 


4 


8 


8 


10 


12 


14 


16 


18 


20 


A H (p) 





-1 


-1 


-4 


-2 


-6 


-3 


-6 


1 


-1 



Table 1. 



However, for p > 22 this difference increases rather rapidly: 



p 


22 


24 


26 


28 


30 


32 


34 


36 


38 


40 


42 


44 


46 


48 


50 


A H (p) 


12 


14 


35 


42 


75 


90 


138 


165 


231 


274 


364 


426 


544 


631 


782 



Table 2. 



Also, an interesting observable phenomenon is a regular oscillation of AJj(p) in 
contrast to the monotonicity of A' c (p). Indeed, in both tables [1] and [2] we have 

signA^(p) = (-If / 2+1 . (6.5) 

for the second difference. This can be conjectured for all p as well as the results of 
observations above. 
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